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ABSTRACT 

A possible discrepancy found in the determination of mass from gravitational lensing 
data, and from X-rays observations, has been largely discussed in the latest years. 
For instance, Miralda-Escude & Babul (1995) have found that the mass estimate 
derived from gravitational lensing can be as much as a factor of 2 — 2.5 larger than the 
mass estimate derived from analysis of the X-rays observations. Another important 
discrepancy related to these data is that X-ray imaging, with some spectral resolution, 
suggest that the mass distribution of the gravitating matter, mostly dark matter, has 
a central cusp, or at least that the dark matter is more centrally condensed than the 
X-ray-emitting gas, and also with respect to the galaxy distribution (Eyles et al. 1991), 
at variance to what is expected from the most accepted models of formation of large 
scale structure. Could these discrepancies be consequence of the standard description 
of the ICM, in which it is assumed hydrostatic equilibrium maintained by thermal 
pressure? In analogy to the interstellar medium of the Galaxy, it is expected a non- 
thermal term of pressure, which contains contributions of magnetic fields, turbulence 
and cosmic rays. We follow the evolution of the ICM, considering a term of magnetic 
pressure, aiming at answering the question whether or not these discrepancies can be 
explained via non-thermal terms of pressure. Our results suggest that the magnetic 
pressure could only affect the dynamics of the ICM on scales as small as ^ 1 kpc. Our 
models are constrained by the observations of large and small scale fields and we are 
successful at reproducing available data, for both Faraday rotation limits and inverse 
Compton limits for the magnetic fields. In our calculations the radius (from the cluster 
center) in which magnetic pressure reaches equipartition is smaller than radii derived 
in previous works. The crucial difference in our models comes from our more realistic 
treatment of the magnetic field geometry, and from the consideration of a sink term 
in the cooling flow which reduces the amplification of the field strength during the 
inflow. In addition the magnetic field calculations arc changed after the cooling flow 
has been formed. 
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1 INTRODUCTION 

Since the work of Loeb & Mao (1994), the possibility of 
explaining the discrepancies on mass determinations, found 
by Miralda-Escude & Babul (1995), via non-thermal pres- 
sure support has been widely discussed (see also Wu & Fang 
1996, 1997; Wu et al. 1998). The discrepancy arises from the 
two most promising techniques to obtain clusters of galaxies 
masses. On one hand, the determination of masses in clus- 
ters of galaxies, via X-ray data, is based on the hypothesis 
that the ICM is in hydrostatic equilibrium with the grav- 
itational potential, using the radial profiles of density and 
temperature. There are uncertainties in the determination of 
temperature profiles, particularly for radii > 1 Mpc, and for 



most systems only a mean emission-weighted X-ray temper- 
ature is available (radial temperature profiles are available 
only for a few clusters ,e.g., Allen & Fabian 1994; Nulsen & 
Bohringer 1995). On the other hand, gravitational lensing 
measures the projected surface density of matter, a method 
which makes no assumptions on the dynamical state of the 
gravitating matter (Fort & Mellier 1994; Miralda-Escude & 
Babul 1995; Small et al. 1997). 

One can find in the literature some attempts to re- 
solve the discrepancy between X-ray and gravitational lens- 
ing mass measurements of clusters of galaxies. For instance, 
Allen (1998) studied in detail a sample of 13 galaxy clusters 
(including cooling fiows, intermediate and non-cooling fiows 
systems) with the goal of comparing X-ray and lensing mass 
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measurements. His conclusions pointed out that, at least for 
cooling flows systems, being more relaxed systems, this dis- 
crepancy is completed resolved, and therefore, non-thermal 
pressures can be discarded in these systems. 

The magnetic field of the ICM can be obtained via Fara- 
day rotation, due to the effect of magnetic field on the po- 
larized radio emission from the cluster or the background 
radio sources. The polarization plane of linearly polarized 
radiation is rotated during the passage through a magne- 
tized plasma. The angle of rotation is <fi = {RM)X^, where 
RM is the rotation measure and A the radiation wavelength 
(Sarazin 1992, for a review). In clusters with diffuse ra- 
dio emission. X-ray observations can give a lower limit to 
the strength of the magnetic field. Typically, this limit is 
B > 0.1 pG (Rephaeh et al. 1987) on scales of ~ 1 Mpc. 
In the case of Faraday rotation the information obtained is 
the upper limit on the intensity of the field, and the mea- 
sured values are RM < 100 rad/m^, that is more or less 
consistent with a intracluster field of B ~ 1 ^G, with a 
coherence length of Ib < 10 kpc. This strength of the mag- 
netic field corresponds to a ratio of magnetic to gas pressure 
oi ps/Pgas < 10~^, implying that B does not influence the 
cluster dynamics (at least on large scales). 

At inner regions, of the cooling flow clusters, the mag- 
netic fields are expected to be amplified due to the gas 
compression (Soker & Sarazin 1990). If they are frozen 
in the mass flow flux, and if this flux is homogeneous 
and spherically symmetric, B ex r~^ and RM oc r^^, 
{Pb oc r^^ and the gas pressure increases slowly). Even 
in this case ps reaches equipartition at a radius tb of 

TB ~ 1 kpc {j^y^ (M/lOO M© yr-l)^''^ In these inner 
regions many sources with very strong Faraday rotations 
were observed, in which the rotation measure can reach val- 
ues of RM ~ 4000 raA/rc? (radio sources associated with 
the central galaxies of the clusters with very strong cooling 
flows (M87/Virgo, Cyg A, Hydra A, 3C 295, A1795)), im- 
plying, B > IQ pG &t Ib ^ I kpc (Taylor & Perley 1993; 
Ge & Owen 1993, 1994). These observations strongly sug- 
gest that the Faraday rotation is created by magnetic fields 
within the cooling flow clusters. 

Another promising method to estimate the cluster 
scale magnetic field, as cited above, is the detection of 
co-spatial inverse Compton X-ray emission with the syn- 
chrotron plasma emission (the 3 K background photons 
scattering off the relativistic electrons can produce a dif- 
fuse X-ray emission) (Rephaeli & Gruber 1988). Therefore, 
this method provides limits on the cluster scale magnetic 
fields, in addition of limits on the non-thermal amount of 
X-ray emission (or even on the relativist electrons energy) in 
galaxy clusters. Such a kind of detection of clusters magnetic 
fields leads, using ROSAT PSPC data and also 327 MHz ra- 
dio map of Abell 85 (a cooling flow cluster, with a central 
dominant cD galaxy and about 100 Me/yr), to an estimate 
of (0.95 ±0.10) pG (Bagchi et al. 1998). However, even non- 
cooling flows clusters present this diffuse, relic radio source 
which can be used to estimate magnetic field strength. For 
instance Ensslin & Biermann (1998) studied limits on the 
Coma cluster magnetic field strength, using this multifre- 
quency observations. They shown that the central magnetic 
field limit is B > 0.3 pG. Others have determined the 
strength of the magnetic field for Coma cluster, using dif- 



ferent techniques and obtaining similar values: B < 1.2 pG 
(Lieu et al. 1996); B > 0.4 pG (Sreekumar et al. 1996). For 
the same cluster (Coma), but using Faraday rotation mea- 
sure, Feretti et al. (1995) estimated magnetic fields of 6.0 pG 
(at scales of 1 kpc), and of 1.7 pG (at scales of 10 kpc) was 
estimated by Kim et al. (1990). 

The above scenario allow us conclude that for both 
methods the observational resolution of the telescope limits 
the detection of smaller scales magnetic fields, implying that 
at scales smaller than 1 kpc the magnetic field strength can 
be higher (Ensslin et al. 1997). Another point to be noted is 
that Faraday rotation measures always gives values higher 
than inverse Compton/CBM measures. Anyway, these fields 
are present in the ICM and therefore justify the such a kind 
of study we present here. Other theoretical works concerning 
the magnetic pressure on the ICM are available (for instance 
Soker & Sarazin 1990; Tribble 1993; Zoabi et al. 1996) and 
we briefly compare our results with those obtained by these 
authors. 

Our goal in this paper is trying to answer the question 
whether or not magnetic support can be relevant in cooling 
flow clusters, using a more realistic treatment of the mag- 
netic field geometric evolution. The scope of the paper is the 
following: in Section 2 we present the hydrodynamical equa- 
tions and the method applied for their solution; Section 3 
describes our models and results compared to the available 
observations; and in Section 4 we discuss our results in the 
light of others obtained in previous works, as well our main 
conclusions. 



2 EVOLUTION OF THE ICM WITH 
MAGNETIC PRESSURE 

The evolution of the intracluster gas is obtained by solving 
the hydrodynamic equations of mass, momentum and energy 
conservation: 
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where «, p, pt, U are the gas velocity, density, total pressure 
and the specific internal energy. The equation of state relates 
U and the temperature. 
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(fcs is the Boltzmann's constant, mn is the hydrogen atom 
mass and p — 0.62 is the mean molecular weight of a fully 
ionized gas with 10% helium by number). The mass dis- 
tribution, M{r), is due to the contribution of the X-rays 
emitting gas plus the cluster coUisionless matter (which is 
the sum of the contributions of galaxies and dark matter 
- the latter being dominant), i.e, M{r) — Mg{r) -f M^iir). 
Mci{r) follows 



pel (r) = po 1 + 
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in which po and a (the cluster core radius) are related to a 
(the hne-of-sight velocity dispersion) via; 9a^ = 47rGa^po- 

The total pressure pt is the sum of thermal and mag- 
netic pressure, e.g., pt = p + Pb- The constraints to the 
magnetic pressure come from observations, from which ps = 
B^/Stt ~ 4 X 10"" erg cm"^ s'^ (cf. Bagchi et al. 1998) for 
a diffuse field located at ~ TOO/Ijq^ kpc from the cluster 
center. Along this paper we will use mostly the ratio be- 
tween magnetic and thermal pressures, or the /3-parameter, 
P = Pb/p- 

The sink term top in the mass equation describes the 
removal of mass from the gas fiow by thermal instabilities. 
The importance of the gas removal was studied in detail 
by Friaga (1993) following the q-description described by 
White & Sarazin (1987). The sink is particularly impor- 
tant when one searches for a steady state solution of the 
cooling flow without an implausible huge accumulation of 
mass at the center. In fact, the condensations formed by the 
sink will probably give rise to stars, planetary bodies or cold 
dense clouds which in turn will constitute a halo surrounding 
the central dominant galaxy. We assume isobaric removal, 
so that the sink does not introduce any additional term in 
the energy equation. Summing up the physics contained in 
this term one can say that the specific mass removal rate is 
Lu = q/tc where the denominator is the instantaneous iso- 
baric cooling time, such that the removal efficiency q relates 
the cooling time to the growth time scale of the thermal in- 
stability in the cooling fiow. We assume q between 1.0 and 
1.5, which are the g- values found to be more consistent with 
the observations (Friaga 1993). 

The cooling function adopted A(T) is the cooling rate 
per unit volume. Since there is no ionization equilibrium for 
temperatures lower than lO*" K, we adopt a non-equilibrium 
cooling function for the gas at T < 10® K (the recombi- 
nation time of important ions is longer than the cooling 
time at these temperatures). The cooling function was cal- 
culated with the atomic database of the photoionization 
code AANGABA (Gruenwald & Viegas 1992). The adopted 
abundances are sub-solar as appropriate for the ICM (Edge 
& Stewart 1991; Fabian 1994; Grevesse & Anders 1989). 

Despite the presence of steep temperature gradients we 
did not consider thermal conduction in our models. This 
can be justified using the fact that on a global scale cooling 
fiow clusters contain cooler gas near the center and hotter 
gas further out. Therefore, the presence of cooling fiows is 
itself a proof that thermal conduction effect is, at least, re- 
duced in the ICM. Models show that thermal conduction 
would erase the observed density and temperature gradi- 
ents in cooling fiows, unless it is inhibited (see, for instance, 
Friaga 1986; David & Bregman 1989). It is well known that 
even weak magnetic field, if it is tangled, can inhibit the 
thermal conduction perpendicular to the field lines. More 
recently it has been argued that electromagnetic instabil- 
ities driven by temperature gradients (or electric currents 
in other situations) also can cause this inhibition in cooling 
fiows (Pistinner et al. 1996), even for non-tangled field lines. 

A spherically symmetric Eulerian code is employed for 
the calculations, which are solved via the finite-difference 
scheme based on Cloutman (1980). The grid points are 
spaced logarithmically, with a grid of 100 cells, with the 
first being 50 pc wide. The innermost cell edge is located 
at 100 pc and the outer boundary at twice the tidal radius 



of the cluster. The artificial viscosity for the treatment of 
the shocks follows the formulation of Tscharnuter & Win- 
kler (1979) based on the Navier-Stokes equation. The outer 
boundary conditions on pressure and density are derived by 
including an outer fictitious cell, the density and pressure in 
which are obtained from extrapolation of power laws over 
the radius fitted to the five outermost real cells. The inner 
boundary conditions are adjusted according to whether in- 
fiow (velocity at the inner boundary is extrapolated from 
the velocities at the innermost cell edges) or outfiow (veloc- 
ity is set zero) prevails locally. The initial conditions for the 
gas are an isothermal atmosphere (To = 10^ K) with 30% 
solar abundances and density distribution following that of 
the cluster dark matter. The evolution is followed until the 
age of 14 Gyr. 

The initial /3 value used here was derived from the mag- 
netic field observations (using, for instance, Bagchi et al. 
1998; Ge & Owen 1993, 1994; Ensslin & Biermann 1998; 
Ruy & Biermann 1998). We assume: frozen-in field; spheri- 
cal symmetry for the fiow and the cluster itself; and that at 
r > re (the cooling radius, see below), the magnetic field is 
isotropic, i.e.. 



_D^ 



Bt/2 = B73 



and Ir — h = I (where Br and Bt are the radial and transver- 
sal components of the magnetic field B and Ir and It are the 
coherence length of the large-scale field in the radial and 
transverse directions). In order to calculate Br and Bt for 
r < Tc we modified the calculation of the magnetic field of 
Soker & Sarazin (1990) by considering an inhomogcneous 
cooling fiow (i.e. Mi ^ M varies with r). Therefore, the two 
components of the field are then given by 
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In our models we take as reference radius the cooling ra- 
dius Tc- In fact we modify the geometry of the field when 
and where the cooling time comes to be less than 10^" yr 
(usually adopted as the condition for the development of 
a cooling fiow). Therefore, our condition to assume a non- 
isotropic field is Uoo = ikBT/2p,mHl^{T)p < 10^° yr. After 
the formation of the cooling fiow, in the inner regions of the 
ICM, the magnetic field geometry is changed, following the 
enhancement of the radial component of the field, due to the 
enhancement of the density. 



3 MODELS AND RESULTS 

In this section we present the results of our models. There 
are four parameters to consider in each one of the models: a, 
the cluster velocity dispersion; po, the initial average mass 
density of the gas; a, the cluster core radius; and /3o, the 
initial magnetic to thermal pressure ratio. We adopted the 
removal efficiency q = 1.5. 

The most important results of our models are shown 
on figures we describe below, for which we assume: a — 
1000 km s~^ and a — 250 kpc. First of all, the evolution we 
follow here is characteristic of cooling fiow clusters and in 
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g cm "^ it rises 
1.5 X 10~^* g cm'^ 



this scenario we discuss the evolution ol the basic thermody- 
namics parameters. Considering the overall characteristics of 
our models, we will compare the results with the very recent 
study based on ROSAT observations of the cores of clusters 
of galaxies, by Peres et al. (1998), focusing on cooling flows 
in a X-rays flux-limited sample (containing the brightest 55 
clusters over the sky in the 2 — 10 keV band). Comparing the 
present models with Peres et al. (1998) deprojection results, 
we see that the central cooling time here adopted as our cool- 
ing flow criterion, e.g. tcoo ^ 10^" yr, is typical for a fraction 
between 70% and 90% of their sample. They also discuss 
briefly the cooling flow age, remembering that in hierarchi- 
cal scenarios for the formation of structures in the Universe, 
clusters are formed by smaller substructures by mergers, and 
therefore the estimation of the cooling flows ages (and the 
cluster ages themselves) is complicated. Anyway they de- 
termine the fraction of cooling flow clusters in their sample 
considering a factor of two in the ages and concluding that 
the fraction do not vary that much (from 13 Gyr to 6 Gyr, 
the fraction varies from 70% to 65%). This allow us conclude 
that our models, which present cooling flows since the cluster 
has the age of ~ 7 — 9 Gyr, are typical for their sample. As 
a matter of fact, the time in which the cooling flow struc- 
ture is formed depends strongly on the initial density we 
adopted. For models with po = 1.25 x 10"^* 
on ~ 9 Gyr, while the models with po 
have it formed on ~ 7 Gyr. We will come back to this point 
later while analyzing the field anisotropy. 

The characteristics of our models are summarized us- 
ing four typical set of initial parameters, and discussing 
some details which came up of the study of a larger 
grid of parameters. Therefore, each model is characterized 
by its position in the (po,/3o) parameter space: model 1 
(po = 1.5 X 10"^** 
1.5 X 10"^** gem 
10 
10-28 gcm-3,/3o = 10-^). 

Figure 1 shows the evolution of density and temperature 
profiles corresponding to model 1, from which the presence of 
the cooling flow on later stages of the ICM evolution and at 
inner regions is remarkable if one notices the steep gradients 
of these quantities. In order to better understand how the 
magnetic field geometry is modified after the cooling flow 
formation, e.g., after the steepness on the temperature and 
density gradients, we follow the evolution of the degree of 
anisotropy, using the concepts previously defined on Section 
2, concerning the geometry of the magnetic field. Hereafter 
we called 'degree of anisotropy' the ratio Bt/Br, noting that 
for the isotropic case it results \/2 and the more anisotropic 
the field geometry the smaller is this ratio. Therefore, we 
present on Figure 2 the evolution of the degree of anisotropy 
since ~ 3.3 Gyr, comparing models 1 and III, in which one 
can see, clearly, that the anisotropy begins decreasing on 
earlier times for models with higher po (model 1) than for 
the ones with lower values of po (model III) . From Figure 2 
we are allowed to conclude that the degree of anisotropy can 
be seen as a sensor of the presence of the cooling flow. In 
another words, the change in the degree of anisotropy can 
be used as another criterion to indicate the epoch, on the 
ICM evolution, in which the cooling flow appear. 

These results can also be discussed in the light of some 
observational works in which the limits to the magnetic field 



,-3 



gem ",/?() ~ 10 ); model II (po — 

"^/3o = 10"^); model 111 (po = 1.25 x 

gcm-^,/3o = 10~2); and model IV (po = 1.25 x 




Figure 1. Evolution of the density and temperature profiles. 
Curves represent early and late stages of the ICM evolution, as 
labeled, for model I. Noting the steep gradients of these quanti- 
ties, i.e., the presence of the cooling flow, at inner regions of the 
cluster, on later stages of the ICM evolution. 



strength on large and small scales of the cooling flow clus- 
ters are given. Following such a kind of observations, as 
previously seen in the introduction section, we chose two 
values of magnetic field strength derived by the authors be- 
low. The first one is presented in Bagchi et al. (1998) who 
estimated, from inverse Compton X-ray emission with the 
synchrotron emission plasma, a cluster-scale (700 kpc) mag- 
netic field strength of (0.95 ± 0.10) pG for Abell 85 (a cool- 
ing flow cluster with a central dominant cD galaxy and M 
~ 100 Mo/yr) . The second one is presented in two papers of 
Ge & Owen (1993, 1994), in which they present and discuss 
rotation measures and the related intensity of the magnetic 
field, giving a range of this intensity at scales of 10 kpc. 
Therefore, our results for the magnetic field strength and 
also for pressures, on large and small scales, are compared 
to the chosen observed ones, in Figures 3 and 4. Reminding 
that the time on which the cooling flow arises is closely re- 
lated to po, one can expect distinct results on the evolution 
of the field intensity from, for instance, model I to model 111. 
However this evolution can be better explained comparing 
model I (Figure 3) to model II (Figure 4), since these two 
models have the same initial density but distinct /3o. 
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Figure 2. Evolution of the anisotropy degree, the tangential to 
radial magnetic field components ratio, for model I (dashed lines) 
and III (full lines), on late stages of the ICM evolution. Noting 
that the presence of the cooling fiow, at inner regions, ^ 200 kpc, 
and on late stages, ^ 9.6 Gyr, of the ICM evolution (see Figure 
1) matches very well with the decrease of the anisotropy degree, 
and that it occurs earlier for model I. 




Figure 3. Evolution of the magnetic strength profiles com- 
pared to the observations. Curves represent early and late stages 
of the ICM evolution, as labeled, for model I (po = 1-5 X 



10 



'^^ g cm '^■ 



/3o 



10 ). Note that the intensity B of the 



magnetic field at 6.7 Gyr is smaller than at 1.6 Gyr, due to the 
fact that at 6.7 Gyr the ICM is in the verge of developing a 
cooling flow, as we can see from the drop in temperature in the 
central region, as shown in Figure 1. In the evolution of the ICM 
before the onset of the cooling flow, the magnetic pressure keeps 
track of the thermal pressure, following the initial conditions for 
Pb/p = /3o < li S'lid the reduction in the thermal pressure just 
after the onset of the cooling flow is reflected in the evolution of 
B. Only after the cooling flow has been established, leading to 
amplification of B, the intensity of the magnetic field will rise to 
high values. 



Our best model, in terms of the magnetic field strength 
compared with observations, is model I (po ~ 1.5 x 
10"^* gcm"^,/3o = 10"^). From Figure 3 it is possible to 
see that on scales of 700 kpc the magnetic field expected 
for the model is higher than the observed one (considering, 
of course, the profile correspondent to redshift zero, or evo- 
lution times on the order of 13 — 14 Gyr), while on scales 
of 10 kpc the model gives a value lower than the observed 
one. Meanwhile, at least on scales of 700 kpc, the situa- 
tion is inverted if one takes a look on Figure 4, for which 
po = 1-5 X 10"^* g cm"^, but /3o = 10"*. Given the uncer- 
tainties characteristics of the observations, we can say that 
our models are in agreement with the magnetic field estima- 
tions available. 

On Figures 5 and 6 we show the magnetic and thermal 
pressures evolution, or in another words, /3-evoIution, for 
models 1 and II respectively, on later times of the ICM evo- 
lution, in order to analyze when and where magnetic pres- 
sure reaches equipartition. Obviously the magnetic pressure 
is compatible with the magnetic field intensities and may be 
compared to the values determined by, for instance, Bagchi 
et al. (1998), ps = B^/Stt ~ 4 x 10"" erg cm" * s"\ at 



scales of 700 kpc, on the present time. From the analysis 
of the magnetic pressures expected from our models it is 
clear that they agree, as well as the magnetic field strength, 
with the observations. Here again model I appears being 
the best one, with (/3o = 10"^), but the values expected 
from model II are not far away from the observed ones as 
well. Noting also that magnetic pressure and/or magnetic 
intensity does not change very much after 12 Gyr, for both 
cases. Results presented on Figures 3-6 would indicate that 
we should adopt an intermediate initial value for /3 (like 
/3o = 5 X 10"'^) in order to obtain a magnetic field intensity 
in better agreement with the observations, at least on scales 
of 700 kpc. Nevertheless such an exercise should not solve 
the match of models and observations on smaller scales, since 
Po — 5 X 10""' should decrease magnetic pressure on scales 
of 10 kpc, at the present time, as a result of the present 
modelling assumptions (see Figure 4). 

Other proposals for the amplification of the magnetic 
field in the center of the cooling flow clusters are: i) ro- 
tational driven mechanisms, in which the twisting of the 
magnetic flux tubes and/or the operation of fast a — cu dy- 
namo are the responsible for the increase of the magnetic 
strength (Godon et al. 1998); ii) turbulence induced am- 
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Figure 4. Evolution of the magnetic strength profiles compared 
to observations. Curves represent early and late stages of the ICM 
evolution, as labeled, for model II (po = 1.5x 10~^® g cm"-^; /3o = 
10~^). See the comments on the evolution of B made in the cap- 
tion of Figure 3. 



plification (Eilek 1990; Mathews & Brighenti 1997). How- 
ever, these processes can not account for the strong magnetic 
fields observed in the center regions, confirming the expec- 
tations previously discussed by authors like Goldshmidt & 
Rephaeh (1993) and Carvalho (1994). 



4 DISCUSSIONS AND CONCLUSIONS 

The present models are in many aspects similar to the one 
of Soker & Sarazin (1990). However there are two important 
differences between our model and theirs: i) they take into 
account only small-scale magnetic field effects; and ii) they 
consider homogenous cooling flow. Since we consider inho- 
mogeneous cooling flow (i.e. M decreases with decreasing r) 
the amplification of B is smaller in our models. As a matter 
of fact the magnetic pressure reaches equipartition only at 
radius as small as <; 1 kpc (model 1) or <; 0.5 kpc (model II), 
because the central increase of the /3 ratio is moderate in our 
model. Our more realistic description of the field geometry is 
crucial. This implies that the effect of the magnetic pressure 
on the total pressure of the intracluster medium, even on 
regions as inner as few kpc, is small. Tribble (1993) study- 
ing the formation of radio haloes in cooling fiow clusters 
from the point of view of the cluster evolution via mergers, 
suggested typical magnetic field strengths of ~ 1 /iG. In ad- 
dition, Zoabi et al. (1996), studying a completely different 
characteristic of the ICM (magnetic fields on the support of 
X-rays clumps and filaments) , adopted the usually assumed 
magnetic to pressure ratio, at few scales of 10 — 20 kpc, of 
0.1, and following a simple geometry of the field in which it 
is amplified by the radial inflow, this ratio become ~ 1 at 
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Figure 5. Evolution of the magnetic (dashed lines) and thermal 
(full lines) pressure profiles on late stages of the ICM evolution 
for model I. Noting that the magnetic pressure increases until 
reaches equipartition at inner regions of the cooling flow (at scales 
^ 1 kpc). 
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Figure 6. Evolution of the magnetic (dashed lines) and thermal 
(full lines) pressure profiles on late stages of the ICM evolution 
for model II. Noting that the magnetic pressure increases until 
reaches equipartition at inner regions of the cooling flow (at scales 
<: 0.5 kpc). 



~ 5 kpc. Again our results are more or less compatible with 
the above ones (for the cluster scale magnetic field), but the 
equipartition condition is reached at smaller scales. 

There are a number of papers discussing heating pro- 
cesses on the inner part of the cooling fiow clusters, in par- 
ticular mechanisms to power the emission lines of optical 
filaments, which use the magnetic energy transformed in op- 
tical emission via magnetic reconnection (Jafelice & Friaga 
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Figure 7. /3 profiles on 14 Gyr for models I - IV. The profiles arc 
quite similars, except for the fact that models with /3o = 10^'* 
have final /3-values lower. From the figure is also clear that the 
equipartition condition occurs at outer radii for higher /3o models, 
and that anyway this condition is reached only on radii smaller 
than ~ 1 kpc. 



1996) or dissipation of Alfven waves (Friaga et al. 1997). 
These works are based in the enhancement of the magnetic 
pressure on scales smaller than ~ 10 kpc, where the fila- 
ments are observed (Heckman et al. 1989). Finally, our re- 
sults suggest that the effect of the magnetic fields on the 
ICM dynamics can be relevant only on very small scales: 
/3 ~ 10~^, r ;$ 10 kpc, and /3 ~ 1, r ;$ 1 kpc, depending on 
the model adopted (see Figure 7). From Figure 7 one can see 
quite clearly that the equipartition condition is reached at 
smaller radii for models in which /3o is equal to 10~^ (model 
II and model IV), emphasizing the agreement between our 
models, another theoretical models, and observations. 

It is also quite relevant noting that the general agree- 
ment of our models and the available data can be empha- 
sized by the fact that observations give us only limits on the 
magnetic field intensities. In the case of rotation measures 
the fimit is the upper one, in contrast with the data coming 
from inverse Compton scattering which give the lower limit 
of this quantity. Therefore, from our best model (model I, 
see Figure 3) the expected field intensity is lower than the 
observed value (provided via rotation measures), on scales of 
10 kpc, and higher than the field intensity derived from X- 
ray inverse Compton scattering, on larger scales (700 kpc). 

That the discrepancy found in the determination of 
mass, from gravitational lensing and from X-rays obser- 
vations (Loeb & Mao, 1994; Loewenstein, 1994; Miralda- 
Escude & Babul, 1995), and in the mass distribution of the 
gravitating matter, mostly dark matter (Eyles et al. 1991), 
can be consequences of the standard description of the ICM, 
in which it is assumed hydrostatic equilibrium driven by 
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thermal pressure (Fabian 1994), is a subject of discussion. 
Allen (1998) argued that, at least for cooling flow clusters, 
the above discrepancy is resolved and, therefore, the effect 
of non-thermal pressures on the hydrostatic equilibrium of 
these systems could be completed discarded. However, it is 
important to point out that the radius in which magnetic 
pressure reaches equipartition is much smaller than the core 
or arc radii obtained by Allen in his analysis (~ 50 kpc, in 
average), implying that, despite Allen's results, at smaller 
scales the non-thermal pressures can be important. 

Theoretical models, like the one here presented, point 
out that magnetic pressure does affect the hydrostatic equi- 
librium of the ICM, but only in the inner radius, as small as 
~ 1 kpc. In addition, it is important to remind that there 
are other sources of non-thermal pressures that could be con- 
sidered jointly to the magnetic pressure before to close the 
discussion on whether or not non-thermal pressures can ex- 
plain the discrepancies on the mass estimations of the galaxy 
clusters. 
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